for ii=1:10000
    ii
check(:,:) = rho(:,:,ii);
[V,D] = eig(check);
dd(:,ii) = diag(D);
vv(:,:,ii) = V;
end